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ABSTRACT 


A monthly mean sea current was calculated for the west coast of the United 
States and the Hawaiian Islands area on a spatial grid of 1° by 1°. These mean 
geostrophic velocities were computer generated from dynamic height data 
obtained from the National Oceanographic Data Center. A method employing 
two-dimensional spline fits of spatially and temporally random hydrographic 


data was developed to determine the monthly averaged geostrophic currents. 


Editor's Note: Reference to a product or comment with respect to it in this 
publication does not indicate, or permit any person to hold out by republica- 
tion in whole or in part or otherwise, that the product has been endorsed, 
authorized, or approved by the U.S. Coast Guard. 
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MONTHLY MEAN SEA CURRENTS FOR THE 
WEST COAST OF THE UNITED STATES AND HAWATI 
BY 
Dr. Kenneth A. Mooney! 
LCDR Alan Summy, MeSo- 
INTRODUCTION 

In 1977 the Oceanographic Unit was given the general task of providing 
improved current information for various areas around the coastal United States 
and Hawaiian Islands. The problem has been that for drift calculations the 
standard method used was to vectorially add a sea current, a wind current, and 
a leeway factor. This method depends on the accuracy of the various components. 
At the time of this study the sea current being used was one derived mainly 
from historic ship set and drift measurements. Set and drift measurements do 
not accurately reflect the sea current both because of the dead reckoning 
navigation errors and the fact that they also contain an inherent wind current 
component. Thus when these quasi sea currents are combined with a wind current 
the result would be to over-drift the object of the search. 

The Oceanographic Unit has endeavored to improve search planning by 
upgrading the mean sea current used. The method used to determine this 
permanent, large scale flow of the oceans required a detailed knowledge of the 
mean temperature, salinity, and density fields in an area. Due to the large 
scope of area to be covered and the prohibitive cost of ship time the 
Oceanographic Unit turned to the National Oceanographic Data Center where such 
information is available for many parts of the world's oceans. 

The result of the analysis of this historical hydrographic data was the 
production of mean sea current charts on a maximum grid of 1° by 1°. These 


mean sea currents are the mean geostrophic currents computed from dynamic 


height anomalies and contain no wind-current contamination. To be used in 
drift forecasting they must be combined vectorially with a wind current 
calculated for the specific time and area of interest. A method for 


calculating a local wind current can be found in Oceanographic Unit Technical 


Report 78-2. 


: U.S. Coast Guard Oceanographic Unit, Building 159-E Navy Yard Annex, 
Washington, D.C. 20593 


PROCEDURES 

The Method 

Since the data was temporally and spatially random it required the use of 
special methods to analyze the data. The most direct method would be to group 
together all the dynamic depth anomaly data for a particular month and then fit 
this data by the method of least squares with a two-dimensional polynomial. 
The virtue of this particular method is that it treats each data point as an 
independent measure of the dynamic depth field. The latitude and longitude are 
used as the independent variable because this is the most convenient form in 
which to obtain the data from NODC. The north-south and east-west velocities 


can now be obtained by solving the following equation: 


U = NO SsreD 

f R sine oo 

C1) 
VS So 1 eb 
fR be 
where 

U = north-south velocity 
V = east-west velocity 
Ff = Coriolis Parameter 
R = Radius of the earth 
@ = Latitude 
dD = Dynamic depth anomaly 
Z~ = Longitude 


The dynamic depth anomaly field frequently displays so much structure that 
to replicate it with a single function would require a polynomial of prohibi- 
tively high order. This difficulty was overcome by the use of two-dimensional 
spline fits of the data. The advantage of spline fitting is that the data 
which was previously fitted by one polynomial is now fitted by several poly- 
nomials, allowing much more structure to be represented. Equation (1) can 
still be used to evaluate the components of geostrophic velocity. 

To put this method into use a spline fit program called SPLAK which 
employs third-order polynomials was obtained from the National Center for 
Atmospheric Research CNCAR) Software Support Library. Even with spline fits, 
the area to be analyzed was too large to be fitted by a single spline fit. Thus 
it was specified that the program have nine rectangles to analyze with the size 
of the rectangles being left as a variable that could be changed as dictated by 
the needs of a particular study area. The reasoning for using nine rectangles 
was that the best fit to the data would occur in the center rectangle where the 
outer boundary effects were least significant. The total area of interest 
would then be analyzed by choosing an appropriate size for the nine rectangles 
according to the amount of structure present and then moving the center rec- 
tangle around the total area of interest with some continuity being supplied 
from one position of the center rectangle to an adjacent position by means of 
the overlapping effect of the eight surrounding rectangles. This approach 
allowed different size areas to be analyzed with equal accuracy and also 
allowed for the development of a generalized program applicable to any area. 
Besides the size of the rectangles, the relative depth to which the currents 
were referenced, i.e. the level of no motion, and the time period over which 
the currents were averaged were both variables which were specified in the 


program. 


For the results presented in this report the program was run with rec-— 
tangles equal to 1° of latitude by 1° of longitude. Since most of the hydro- 
graphic data extended only to 500 meters that depth was chosen as the reference 
level enabling the use of most of the available data. While 500 meters may 
appear to be somewhat shallow for those who normally use 1000 meters or 1200 
meters as a reference level, cross-checking when possible showed no significant 
differences in surface velocities when 500 meters versus 1000 meters was used 
as a reference level. The spline fit was evaluated at the center of rectanges 
%° by 4°, 3° by 3°, and 1° by 1° on their sides. In the final analysis the 
current structure was such that no areas required rectangles smaller than 3° by 
3° to properly characterize the currents. 

Data Presentation 

The results of the analysis are depicted for each area in twelve 
chartlets, one for each month. At the bottom of each block you will note that 
on the left of the slash is the current direction in degrees true and to the 


right of the slash is the current speed in knots. 
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